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SUMMARY 


A method is described for calibrating airspeed installations on 
airplanes at transonic and supersonic speeds in vertical— plane maneuvers 
in which use is made of measurements of normal and longitudinal acceler- 
ations and attitude angle. The method involves starting or ending a 
calibration run near level flight at a speed for which the airspeed 
calibration is known and hence for which the free— stream static pressure 
may be determined. Integration of vertical velocity, computed from the 
normal and longitudinal accelerations and the attitude angle, determines 
the change in altitude which, when combined with the temperature 
measurements, gives the change in free— stream static pressure from the 
start or end of the calibration run and hence the variation of free— 
stream static pressure during the calibration run. In this method, all 
the required instrumentation is carried within the airplane. 

A study of the effect of the various sources of error on the 
accuracy of a calibration indicated that the quantities can be measured 
accurately enough to insure a satisfactory calibration. 


INTRODUCTION 

A method of calibrating the static— pressure source of a pitot— static 
installation on an airplane at high speed and high altitude by the use of 
radar— phototheodolite tracking equipment was described in reference 1. 

In this method, the radar— phototheodolite equipment is used to establish 
the geometric altitude of the airplane in surveys of atmospheric pressure 
made at speeds for which the airspeed calibration is known and in 
maneuvers under conditions for which the calibration is desired. Although 
this method is precise, its use is limited since radar— phototheodolite 
equipment is not generally available. In order that airspeed calibrations 
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may be performed universally, tvo additional methods were devised at 
the Langley Aeronautical Laboratory. One method, described in detail in 
reference 2, makes use of measurements of total pressure, static pressure, 
and temperature. In this method; a survey is made of static pressure and 
temperature at airplane speeds at which the airspeed calibration is known. 
The airplane is then flown through the region surveyed under conditions 
for which the calibration is desired and measurements of total pressure, 
static pressure, and temperature are repeated. The values of total 
pressure and temperature at a given instant in the calibration run, 
together with the calibration constant of the thermometer, are used to 
determine the relation between ambient temperature and static pressure. 
This relation, when compared with that for the survey, determines the 
free— stream static pressure at that instant. This method requires 
precise measurements of temperature. Another method, described herein, 
makes use of measurements of acceleration and attitude angle together 
with measurements required for the previous method. The accuracy of 
this method is determined principally by the accuracy of the measurements 
of normal and longitudinal accelerations and of the attitude angle.. 
Measurements of temperature need not be so precise as those for the 
previous method. The method is restricted to vertical— plane maneuvers. 

SYMBOLS 


P 

Pm 


P- 


A Pm 


t 

^m 


e t 


free-stream static or atmospheric pressure 
indicated static pressure 

density 

(V m ~ 

static— pressure error I — — 

\ P 

change of static pressure indicated by statoscope recorder 

free— stream total pressure for subsonic flow and total pressure 
behind normal shock for supersonic flow 

elapsed time 

measured elapsed time. 

error in elapsed time 
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h altitude 

M free— stream Mach number 

M’ indicated free— stream Mach number 

T free— stream temperature, absolute units 

T m measured temperature, absolute units 

T* temperature defined by equation (5) 

e T error in free-stream temperature 

K temperature recovery factor 

a x longitudinal acceleration, positive forward along x— axis 

a z normal acceleration, positive upward along z— axis 

vertical acceleration, positive upward along vertical of earth 
(defined by equation (10)) 

a^. lateral acceleration, positive to right of y— axis 

a angle between rays of sun and longitudinal axis as measured by 

sun camera 

P longitude of airplane 

7 ratio of specific heats (1.4) 

e declination of sun 

X elevation angle of sun 

6 attitude angle, positive below horizon 

t latitude of airplane 

0 angle of bank, positive right wing down 

i|r angle of yaw, positive to right 


Tm~T \ 

,0.2M 2 T/ 
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to Greenwich hour angle 

v vertical velocity 

R gas constant 

So error in estimating latitude of airplane (distance along 

longitude (3) • 1 

s T error in estimating longitude of airplane (distance along 

latitude t) 

g acceleration due to gravity 

A prefix denoting error 

Subscript: 

1 near start or end of maneuver 


METHOD 


The principal feature of most procedures for calibration of air- 
speed installations on airplanes is the determination of ambient or 
free-stream static pressure. In the method described herein, free- 
stream static pressure is obtained over a range of altitude with the aid 
of the relation that the rate of change of static pressure with altitude 
is equal to the density, or 


dp 

— = -gp 
dh 


( 1 ) 


Since 




dp _ p 

• dh ~RT 


then 
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or 


dp = - JL dh (2) 

RT 

Equation (2) may be integrated as 




(3) 


Therefore, if the pressure at one altitude is known, the pressure at 
other altitudes may be determined provided the ambient temperature and 
the change in altitude are determined. The ambient temperature T may 
be determined from the measured temperature T m with the use of the 

relation 



(V) 


Since only the indicated Mach number M* is known for flight conditions 
where no airspeed calibration exists, the ambient temperature is deter T 
mined only approximately as ■ 

t* = 3a . (5) 

1 + 7 - — 1 KM» 2 
2 


The use of T* in equation (3) would result in a small error in p and 
hence two or more approximations may be necessary. 

An alternate integral form of equation (2) is 



= 1 


— n 



( 6 ) 
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After substitution of the measured pressure and the temperature T* 

in the right side of equation (6), the equation becomes 



The value of n may be selected for the flight conditions encountered 
so that only one approximation is required in the determination of p. 
(See appendix A.) For values of temperature recovery factor K of the 
thermometer near unity and for subsonic and low— supersonic airspeeds, a 

value of n of ~ — or 0.286 gives satisfactory results . 

7 

The change in altitude dh in equation (7) may be determined from 
vertical velocity computed from measurements of pressure and temperature 
at an instant in the calibration run when the airspeed calibration or . 
the static pressure is known and from vertical acceleration computed 
from measurements of longitudinal and normal acceleration and attitude 
angle, as 



and 

a v = a z cos 0 - a x sin 0 — g ( 10) 
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The method as outlined here is limited to maneuvers in vertical planes. 
The method is further illustrated hy a hypothetical calibration in 
appendix D. 

The relation given by equation (7) may also be used to advantage 
in the radar method of reference 1, particularly for tests which would 
not permit a survey of static pressure to be made over the desired range 
of altitude at flight conditions (Mach number and lift coefficient) at 
which the airspeed calibration was known. For such a flight condition, 
however, the airspeed calibration at least at one instant during the 
test must be known and the change in altitude from the altitude at that 
instant would, of course, be determined from radar measurements. 

EQUIPMENT 


The airplane on which the pitot— static installation is to be 
calibrated should be equipped with a recording altimeter, an airspeed 
recorder, a statoscope with a sensitive pressure recorder, a thermometer 
and temperature recorder, an accelerometer with normal and longitudinal 
components, an attitude recorder, and a timer. The recording altimeter 
is used to record the static pressure measured by the static— pressure 
source. The airspeed recorder is used to record the impact pressure, 
or the difference between total pressure and static pressure, measured 
by the pitot— static installation. The statoscope with the sensitive 
pressure recorder and timer is used to measure accurately the time rate 
of change of static pressure for the determination of vertical velocity 
near the start or end of the calibration run. The airspeed recorder, 
the recording altimeter, and the statoscope should be the only instruments 
connected to the static-pressure source and should be located as near as 
possible to it in order to minimize the lag of the pressure system. 
Provision must also be made so that the large volume of the statoscope 
is not continuously open to the static— pressure source because of lag 
considerations and that the sensitive pressure cell is not subjected to 
pressures beyond its limit in order to prevent damage. The magnitude 
of the pressure lag of the airspeed installation may be determined by 
methods described in reference 3 . Where the lag is appreciable, 
corrections must be made to the measured static pressure. The thermometer 
is used to determine free— stream temperature'. A properly shielded 
thermometer with a high temperature recovery factor (approaching 1 . 0 ) 
is recommended since it should be least affected by position on the 
airplane. The thermometer should also have a low lag; corrections 
should be made if the lag is appreciable. The attitude recorder is used 
to measure the attitude of the airplane and may consist of a horizon 
camera, a sun camera, or an attitude gyroscope. A horizon camera 
shooting either forward or laterally is probably the most desirable 
attitude recorder in localities where the horizon is not obscured by 
haze. The attitude gyroscope and the sun camera, however, may be more 
generally used. The attitude gyroscope measures the change in attitude 
angle. The attitude angle at some instant during the calibration must 
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therefore he determined from other measurements, perhaps from the 
statoscope measurements and from the estimated angle of attack. When 
a sun camera is used, the airplane should be equipped with an indicating 
device to enable the pilot the fly the airplane in a vertical plane with 
the lateral axis normal to the rays of the sun. One such device is a 
sun dial. The determination of the attitude angle from measurements 
with the sun camera is discussed in appendix B. 


CALIBRATION PROCEDURE 


The calibration procedure described herein may be used for level 
flight, climb, dive, push-down, and pull-out in a vertical plane. A 
calibration may be obtained in a single run for some conditions and does 
not require any additional survey of static pressure and temperature as 
do the methods of references 1 and 2. 

For a dive starting with a push-down from level flight, the airplane 
is flown prior to the push-down at a condition (lift coefficient and Mach 
number) for which the airspeed calibration is known. Also, prior to the 
push— down, the statoscope is closed and the recording instruments are 
turned on. Measurements are made throughout the calibration run of 
impact pressure, static pressure, temperature, longitudinal acceleration, 
normal acceleration, attitude angle, and the change in static pressure 
from the pressure at which the statoscope was closed. The differential- 
pressure recorder of the statoscope, however, may go off scale when the 
differential pressure exceeds the range of the sensitive pressure cell. 
When the dive includes a pull-out to level flight, or a pull-out, 
climb, and then a push- down to level flight, it may be desirable, but not 
necessary, to obtain further measurements with the statoscope. The pull- 
out or push-down is made to a speed for which the airspeed calibration 
is known. The statoscope is then opened momentarily and closed. The 
instruments are turned off a few seconds after the statoscope is closed. 

For a dive starting with a half roll and a pull-out (split— S) the 
instruments are turned on after the airplane has attained the desired 
attitude in the- vertical plane and are kept on throughout the pull-out 
to level flight, or throughout the pull-out, climb, and push— down to 
level flight, whichever the case may be. The pull-out or push— down is 
made to a speed for which the calibration is known and the statoscope 
is opened momentarily and then closed. The instruments are turned off 
a few seconds after the statoscope is closed. 

For a climb starting from level flight, the procedure is essentially 
the same as that for the dive starting from the push-down. 
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When a sun camera is used to obtain the attitude angle, the airplane 
must be flown with the lateral axis normal or very nearly normal to the 
rays of the sun. A simple sun dial mounted ahead of the canopy may be 
used by the pilot as an indicator for keeping the lateral axis normal to 
the rays of the sun. The local time when the instruments are turned on 
should be determined accurately. .The. clock that is used for this purpose 
may be checked against radio time signal (National Bureau of Standards 
radio station WWV). 

When a horizon camera or an attitude gyroscope is used, the 
airplane must be flown in a vertical plane. The airplane, however, is 
,not restricted to any particular vertical plane as in the case with the 
sun camera. 


ACCURACY 


The accuracy of the method depends principally on the accuracy of 
determining the attitude angle and the accuracy of the accelerometer 
measurements. The equations for the various errors are derived in 
appendix C and the errors in the computed quantities due to assumed 
errors in measured quantities are presented in figures 5 to 15 . When 
a sun camera is used to determine the attitude angle, errors in attitude 
angle may arise from an error in time of the calibration, an error in 
longitude and latitude of the airplane, and an error in measurement of 
angle with the sun camera. By measuring the time of the start or end 
of the calibration run to within a few seconds, the error in attitude 
angle due to an error in the measurement of the time of the calibration 
may be almost eliminated (see fig. 5). The clock should be checked 
against an accurate source of time, perhaps against radio t ime signal, 
before or after the calibration. Because of the speed of the airplane, 
and hence the distance covered in a calibration run, the latitude and 
longitude of the airplane may vary appreciably during the run. If, 
however, the pilot can estimate his position to within 10 miles, the 
error in attitude angle can be kept to within ±0.1° (fig. 6). The 
error in the angle of sun rays relative to the longitudinal axis as 
measured with the sun camera need not be more than about 0.1° for a 
properly designed sun camera. The probable maximum error in attitude 
angle should be of the order of 0.1°. 

- A properly designed horizon camera can probably measure attitude 
angle to within 0.1°. When the horizon camera is shooting forward, the 
correction for dip of the horizon can be estimated to within +0.1°. 

When the attitude angle 'is obtained from the flight— path angle and 
the angle of attack at one instant in the calibration run, together with 
the change in attitude angle as measured with an attitude gyroscope for 
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other times during the calibration run, the error in attitude angle 
may be of a larger magnitude than that for the method with which the 
sun camera is used. The error in the flight-path angle arises from an 
error in determining vertical velocity with the use of the statoscope. 

If the rate of change of static pressure with time is determined to an 
accuracy of ±0.01 inch of water per second, the error in vertical 
velocity at ko, 000-feet altitude is about 3 feet per second. This error 
in vertical velocity for a flight Mach number of about 0.6 corresponds 
-to an error of about 0.3° in flight— path angle. The error in esti m ating 
the angle of attack at the instant when the flight— path angle is deter- 
mined depends on the applicability and the accuracy of the information 
on which the estimate is based. The NACA attitude gyroscopes have a 
drift, of about 3° P er minute. For a calibration run lasting, for 
example, 20 seconds, the error in attitude angle due to the drift of 
the gyroscope would be of the order of 1° at the end of the calibration 
run. The probable maximum error in attitude angle (assuming a 
0.3° error in estimating the angle of attack) for a calibration run 
lasting 20 seconds would vary from about 0.4° at the start to about 1.1° 
at the end. 

The error in the vertical component of acceleration due to an error 
of ±1° in attitude angle as shown in figure 7 increases with increase in 
normal acceleration and with increase in attitude angle. In dives, 
normal acceleration would vary from about l.Og near level flight to 0 in 
a vertical dive. The error in vertical acceleration in a dive would 
probably be of the order of O.Olg for an error of ±1° in attitude angle. 
In a pull-out, the error in vertical acceleration would be larger but 
would probably occur only near the end of the calibration run. 

The error in the vertical component of acceleration due to an error 
of ±0.01g in the normal and longitudinal components of acceleration is 
of the order of O.Olg (fig. 8). An NACA recording accelerometer has an 
accuracy of l/b percent of full range. A longitudinal accelerometer of 

i|g range and a normal accelerometer of 0 to 1 g range would have an 

accuracy of 0.0025g. A normal accelerometer of 0 to kg range would have 
an accuracy of O.Olg. In order to get improved accuracy for normal 
accelerations, a combination low— range and high— range normal acceler- 
ometer can be used. For example, in a calibration run involving a dive 
and pull-out, a low— range normal accelerometer can be evaluated for the 
dive and the high-range accelerometer for the pull-out. 

The error in the vertical component of acceleration due to 
neglecting the angle of bank varies with attitude angle and normal 
acceleration (fig. 9). For. l.Og normal acceleration, neglecting an 
angle of bank of 10° results in an error in vertical acceleration of 
-0.015g at zero attitude angle and of Og in a vertical dive. 
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The error in the vertical component of. acceleration due to 
neglecting the angle of yaw is shown in figure 10. For a 2° angle of 
yaw, the error in vertical acceleration is less than O.Olg. For transonic 
and supersonic speeds, an angle of yaw of 2° probably would not be 
unintentionally exceeded and certainly would not be maintained. 

The error in static pressure due to a constant error of ±0.01g in 
vertical acceleration varies as the square of the time for the cali- 
bration run (fig. 11). For a calibration run lasting 20 seconds, the 
error, as percent of free— stream static pressure, is O.3., and for 
40 seconds, the error is 1.1. 

The error in vertical velocity due to an error in time rate of 
change of static pressure as determined from the statoscope is shown in 
figure 12. For an altitude of 40,000 feet, the error in vertical 
velocity is 2.8 feet per second for an error of ±0.01 inch of water per 
second in the time rate of change of static pressure. 

The error in static pressure due to an error of ±0.01 inch of 
water per second in the time rate of change of static pressure is given 
in figure 13 as a function of time for an altitude of 46,000 feet. The 
errors in static pressure given as percent of static pressure are about 
0.3 and 0.5 for calibration runs lasting 20 and 40 seconds, respectively. 

f s 

The error in static pressure due to an error of ±1 percent in 
measured temperature (or about ±5° F) varies with the range of pressures 
covered in the calibration (fig. l4) . For a range of static pressures 
from 0.7 to 3 times the initial static pressure, the error is within 
1 percent of free— stream static pressure. 

An error in total pressure results in an error in Mach number 
and, hence, in static pressure as determined from equation (7)* 

Errors in total pressure due, for instance, to angularity of flow may 
be avoided by the use of properly designed pitot tubes . The error in 
free— stream static pressure due 1 to an error of ±1 percent of total 
pressure is seen in figure Ik to be within 0.3 percent of free— stream 
static pressure for a range of static pressures from 0.7 to 3 times the 
initial static pressure. 

In evaluating the free— stream static pressure with the use of 
equation (6), the value of the gas constant may be taken to be 53- 3y 
the value for dry air. The resulting error due to neglecting moisture 
content may be shown to be negligible. For example, if the air were 
saturated, the use of the gas constant for dry air would introduce an 
error of about 0.2 percent of free— stream static pressure in a dive 
from an altitude of 10,000 feet to sea level and of 0.01 percent in a 
dive from 30,000 feet to 20,000 feet. 
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The error in static pressure due to error in an NACA recording 

altimeter is ±i percent of the full-scale reading. For an altimeter 

covering a range from sea level to 50,000 feet, the error in static 
pressure is ±1 inch of water or ±1.3 percent of static pressure at an 
altitude of 40,000 feet. For an altimeter covering a range of altitudes 
above 30,000 feet, the error is ±0.2 inch of water or ±0.3 percent of 
static pressure at 40,000 feet. Further improvement in accuracy of the 
static-pressure measurements may be obtained with the use of a statoscope 
equipped with a differential— pressure recorder having a range to cover 
the change in static pressure over the desired range of altitudes. The 
error in Mach number due to an error of ±1 percent of static pressure 
is shown in figure 15 . 

The elapsed time t may be measured to within 0.01 percent with 
the, use of a tuning-fork timer. The static-pressure error corresponding 
to this error in time is (according to equation (^5)) within 0.05 per- 
cent of free-stream static pressure for a range of altitudes (h — hq) 
of 50,000 feet. Since all instruments can be of the continuous recording 
type, no consistent error should result from the correlation of these 
records . 

When a calibration run begins and ends near level flight at a 
speed for which the calibration is known, a check on the constant 
errors in calibration (due to errors in acceleration or attitude angle) 
is obtained in that the vertical velocity at the end of the calibration 
run should be equal to vertical velocity at the start plus the time 
integral of the vertical acceleration. A consistent error in the 
recording altimeter does not affect the determination of the static- 
pressure error since the consistent error would be included in the 
computed free-stream static pressure as well as in the measured 
pressure . 

For the hypothetical airspeed calibration illustrated in 
appendix D, the probable maximum error in the determined static— pressure 
source was estimated on the basis of the following assumed accuracies: 

±3.2 
±0.0025 
± 0.01 

± 1.0 

±3.01 


A 0, degrees 

Aa i (for -0.5g to 0.5g longitudinal accelerometer), g units . 
Aa n (for-0 to 4g normal accelerometer), g units 

Ap ( for altimeter with a range including sea— level pressure) , 
inches of water 

/dp \ 

A — , inches of water per second - 

AT m , °F 


± 1.0 
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The probable maximum error in the determined static— pressure error of 
the static— pressure source at the end of the dive was estimated to be 
+1.2 percent of free— stream static pressure. If a normal accelerometer 
having 0 to 1 g range and an altimeter measuring pressures only above 
an altitude of 25,000 feet are used, the probable maximum error is 
reduced to +0.6 percent of free—stream static pressure. 


CONCLUDING REMARKS 


A method is described for calibrating airspeed installations on 
airplanes at transonic and supersonic speeds in vertical— plane maneuvers 
in which use is made of measurements of normal and longitudinal acceler- 
ations and attitude angle. The method involves starting or ending a 
calibration run near level flight at a speed for which the airspeed 
calibration is known and hence for which the free— stream static pressure 
may be determined. Integration of the vertical- acceleration computed 
from the normal and longitudinal accelerations and the attitude angle 
determines the change in altitude which, when combined with the temper- 
ature measurements, gives the change in static pressure from the start 
or end of the calibration run and hence the variation of free-^stream 
static pressure during the calibration run. The static— pressure error 
is then obtained at any instant during the calibration run by subtracting 
the free— stream static pressure from the indicated static pressure. 

In the method described herein the required instrumentation is 
carried within the airplane . Should the airplane at any time enter a 
previously unexplored flight condition in a vertical— plane maneuver, a 
calibration may be readily obtained. 

In measuring the attitude angle, the sun camera or the horizon 
camera appears to offer better accuracy than the use of an attitude 
gyroscope. The airplane with the sun camera, however, would be limited 
to a vertical plane in which the lateral axis of the airplane is normal 
to the rays of the sun, whereas the airplane with either the horizon 
camera or the attitude gyroscope could select any vertical plane. The 
horizon camera is restricted to use In localities where the horizon is 
not obscured by haze . 
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A study of the effect of various sources of error on the accuracy 
of a calibration indicated that the required quantities can be measured 
accurately enough to insure a satisfactory calibration. 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va . , February 23, 1950 
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APPENDIX A 

CALCULATION OF SUITABLE VALUES OF n FOR USE IN EQUATION (7) 


The value of n that yields zero or nearly zero error in the 
computed free— stream static pressure as a result of using Pju and M' 

in equation (7) may be found by first differentiating equation (6) as 



( 11 ) 


where e p in the integral is the error in the static— pressure source 

and Ap/p on the left side of the equation is the error in the computed 
free— stream static pressure due to use of p m and M f . Also in the 

integral e T is the error due to use of M* . For zero error in 
computing free— stream static pressure 

ne p = e T 


or 


n = 


For M 1.0, the value of — may be determined from equation (4) 

e p 

and the equation 


PT = p (i * 2_i 1 


( 12 ) 
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Differentiating each equation and combining results in the following 
expression: 


n = 


= p dT 
e p T dp 


7 






K 


(13) 


For M 1.0 the value of 
equation (4) and the equation 


n may be similarly computed from 


7 + 1 „2 
PT = — o — M P 


(7 + 1)^M‘ 


2 m 2 


47M 2 — 2(7 — l) 


1 

7-1 


(Ik) 


or 


n 


7-1 27 M 2 - (7 - 1) M 2 K 


27 


2M 2 - 1 q + km 2 
2 


(15) 


The values of n for K = 0.9 and K = 1.0 are tabulated for various 
Mach numbers as follows : 


M 

n 

11 

0 

VO 

K = 1.0 

0.5 

0.258 

0.286 

1.0 

.261 

.286 

1.5 

.347 

-=t 

r- 

m 

. 2.0 

.461 

.490 

3.0 

.644 

.670 


.759 

. ..... .... . j 

.779 
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Since, in the general case, a range of Mach number would be covered in 
a calibration test, a mean value of n may be taken. The possible 
error in the computed free— stream static pressure resulting from the 
use of the mean value for a range of Mach numbers from M = 1.0 
to M = 2.0 was estimated to be less than 3 percent of the static- 
pressure error of the airspeed installation. 
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APPENDIX B 

CALCULATION OF ATTITUDE ANGLE FROM MEASUREMENTS WITH SUN CAMERA 


The attitude angle is determined by subtracting the elevation 
angle of the sun from the angle recorded by the sun camera 


9 = a — X 


(16) 


The elevation angle of the sun may be found in navigational tables by 
use of the date of calibration, the time at the start or end of the 
calibration rim, and the longitude and latitude of the airplane. The 
elevation angle may also be found from the expression 


sin X = sin € sin r + cos e cos t cos(<d — p) 


(17) 
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APPENDIX C 

CALCULATION OF ERRORS 


Error in Attitude Angle Due to Error in Time of Calibration 

An error in the time of calibration results in an error in the 
elevation angle of the sun and hence in the attitude angle. The error* 
in the elevation angle of the sun is, after differentiation of 
equation (17), 


cos c cos t . f . 

AA = sm(o) — p)A03 

cos A 


(18) 


or 


COS 6 COS T . / 15 

AX sin(a> — 3) — — At 

cos X 3600 


(19) 


where 

Ato = At 

3600 

and At is the error in time in seconds . 
Since 


A X — — r\G 


A 6 = 0.004l6 — g C0S T sin(o) — 3) At (20) 

cos X - 

The error in attitude angle due to an error of 60 seconds in time of 
the calibration is plotted in figure 5- 
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Error in Attitude Angle Due to Error in Latitude 
and Longitude of the Airplane 

The error in attitude angle due to error in the latitude of the 
airplane is 


AQ = — AX — 


—sin e cos t + cos e sin t cos(o> — p) 


cos X 


At (21) 


or 


A 9 = 


-sin e cos t + cos e sin t cos(ab — g) 97.3 
cos X J i}000 


( 22 ) 


where 


at 


4000 13 


and Sp is in miles. 

Similarly, the error in attitude angle due to error in the 
longitude of the airplane is 


. A0 = -<os € sin(co-p) 2L3 St (23) 

cos X 4000 


where St is in miles. 

The error in attitude angle due to error in the latitude and 
longitude of the airplane is plotted in figure 6. 


Error in Vertical Component of Acceleration 
Due to Error in Attitude Angle 

The error in vertical acceleration due to an error in attitude 
angle is, after differentiation of equation- ( 10) , 
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ABry = — (a 2 sin 0 + cos 0)A0 


(24) 


where A0 is in radians. The error in the vertical component of 
acceleration for a ±1° error in attitude angle is shown in figure 7* 


Error in Vertical Component of Acceleration Due to Error in 

Longitudinal and Normal Components of Acceleration 

The error in the vertical component of acceleration due to error 
in the longitudinal component of acceleration is, after differentiation 
of equation (10), 


. A&y = -i sin 0 


(25) 


Similarly, for an error in the normal component of acceleration 


Ahy = Aa z cos 0 


( 26 ) 


The error in the vertical component of acceleration due to an error 
of ±0.01g in the longitudinal and normal components of acceleration is 
presented in figure 8. 


Error in Vertical Component of Acceleration Due to Angle of Bank 

The expression for vertical acceleration, including the angle of 
bank, is 

a^- = a z cos 0 cos 0 — a x sin 0 — g - (27) 

where 0 is, the angle of bank. 

If the angle of bank is neglected, the error is 

= a z cos 0(cos 0 — 1) ( 28 ) 
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The error in the vertical component of acceleration is plotted in 
figure 9 for angles of bank of 5°, 10°, and 15°. 


Error in Vertical Component of Acceleration Due to Angle of Yaw 

The vertical component of acceleration, including the effect of 
yaw, is 


a v = a z cos 9 — a x sin 9 cos + a y sin i|r sin 9 — g (29) 

where i is the angle of yaw and a y is the lateral acceleration due 
to yaw. 

If the angle of yaw is neglected, the error is 

Aay = a x sin 0(l — cos i) + a y sin sin 9 

= ja x (l — cos i|f) + ay sin tjsin 9 (30) 

This error in the vertical component of acceleration is shown in 
figure 10 for 2° and 4° of yaw. 


Error in Free-Stream Static Pressure Due to Error in 
Determining Vertical Component of Acceleration 

The error in free— stream static pressure due to error in the 
vertical component of acceleration is found by substituting e ques- 
tion (8) into equation (7) and then differentiating the resulting 
equation, or 


£p = _1_ 

P "pn 



p n (l + 0.2KM2) 


RT, 


m 


dt 


dt 


(31) 


In the integral My, p, M, and T m may vary with time ^ In order to 

obtain the order of magnitude of the error, however, it is sufficient 
to assume these quantities as constant. Then 



NACA TN 2099 


23 


AP _ _^v +2 

p ” 2RT 


( 32 ) 


The variation of the error in static pressure with time is shown in 
figure 11 for an error in vertical acceleratipn of +0.0ig. 


Error in Vertical Velocity Due to Error in Determining 

the Time Rate of Change of Static Pressure 

The error in. vertical velocity due to error in determining the 
time rate of change of static pressure is obtained from equation (9) 
as 




(33) 


The error in vertical velocity due to an error of 0.01 inch of water 
per second in the time rate of change of static pressure is shown in 
figure 12 for various altitudes. 


Error in Free-Stream Static Pressure Due to Error in 

Determining Vertical Velocity or the Time 

Rate of Change of Static Pressure 

The error in free— stream static pressure due to error in 
determining vertical velocity is, after substitution of equation (8) 
into equation (7) and differentiation of- the resulting equation, 


_ __1_ P 1 p n (l + 0.2KM 2 ) 
p P 11 Jt! RT m 


(3b) 
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In evaluating the order of magnitude of error in static pressure, T^, M, 
and p may he assumed as constant in the integral and therefore 


. Ap _ Avt 

■p" “ “RT~ 


(35) 


In terms of error in time rate of change of static pressure, the error 
in static pressure is, after substitution for Av from equation (33)> 


An _ a / dp\ t 
p \dt/p 


(36) 


This static— pressiire error is plotted in figure 13 as a function of time 
for an altitude of about 40,000 feet and an error of ±0.01 inch of water 
in the time rate of change of static pressure. 


Error in Free-Stream Static Pressure Due to Error .in Measuring T m 

The error in free— stream static pressure due to error in 
measuring T m is, after differentiation of equation (7), 


£p _ 1 P 1 * p n (l + 0.2KM 2 ) AT m 


P U , PI 


RTi 


m 


T m 


dh 


(37) 


AT m 

For = Constant, equation (37) reduces to 

T m 


Ap _ A^m 
P nT, 


m 


Pi 

— I - 1 


(38) 


The error in free— stream static pressure due to an error of ±1 percent 
in T m is plotted in figure 14 against the ratio of initial pres- 
sure to pressure p at any time during the calibration.. A value of 

n of or 0.286 was assumed. 


7 
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Error in Free— Stream Static Pressure Due to 
Error in Measuring Total Pressure 

The error in total pressure introduces an error in the computation 
of temperature and hence in the computation of free— stream static pres- 
sure. The static-pressure error may he found by differentiating 
equation (6) as 


ap = 1 / h at 

P " P* A-l RT T 


(39) 


AT Ap T 

If — is assumed to be a constant and is related to through the 

T • • P T 

use 'of equations (4), (12), and (l4), the static— pressure error is 
found to be 


Ap _ apt /piY 1 _ 1 

p =_p t L • 


m 


■The error in static pressure due to ±1— percent error in total pressure 
is plotted in figure l4 against the ratio of initial pressure p^ to 

pressure p at any time during the calibration for n = 2 — ■ . 

7 


Error in Free— Stream Static Pressure Due 
to Error in Elapsed Time 

Integration of equation (8) between the limits of h and h^ 
yields 



w 


t 
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If the measured elapsed time tm is substituted into equation (4l) the 
equation becomes 


K ~ h l = 





If the error in elapsed time is defined as 


e t 


t 


the error in altitude determined from equations (4l) and (42) becomes 



(43) 


The corresponding error in free— stream static pressure is 


= fi 
p RT 


v^t — 2(h — h^) 


(44) 


The error in free— stream static pressure , is a maximum if v-^ is zero 
or if v^_ has a direction opposite to the resultant change in altitude. 
For v-^ = 0, 


Ap 2e t ( h x - h) 

T = RT 


m 


Error in Mach Number Due to Error in Free— Stream Static Pressure 


The error in Mach number due to an error in static pressure is, 
after differentiation of equation (12) for M ^ 1.0, 


AM = - 


1 + 0.2M 2 

n 


p 


(46) 
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and, after differentiation of equation (l4) for M ^ 1.0, 


AM = 


M 7M 2 - 1 AP 
7 2M 2 - 1 P 


( 47 ) 


The error in Mach number due to a ±1— percent error in free— stream static 
pressure is shown in figure 15* 
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APPENDIX D 

ILLUSTRATION OF METHOD WITH THE USE OF HYPOTHETICAL DATA 


In the hypothetical airspeed calibration, the airplane is assumed 
to be equipped with a sun camera for the determination of the attitude 
angle. The following data, pertinent to the determination of the 
elevation or altitude angle of the sun, are assumed: 

Date April 25, 19^9 

Time at start of calibration lS^OO 111 E.S.T. or lS^OO 111 G.C.T. 

Longitude of Langley Laboratory, (3 76°21' 

Latitude of Langley Laboratory, t 37°5* 

From navigational tables, the Greenwich hour angle, co is found to 
be 90°31* and the declination e is 13°15' . With these values of 
longitude, latitude, Greenwich hour angle, and declination, the 
elevation angle of the sun A is found from navigational tables or 
computed from the expression (see appendix B) 


sin A = sin'e sin t + cos e cos t cos(o> — (3) 
to be 65. 4°. The attitude angle of the airplane is then obtained as 


Q = a — A 


= a — 65. 


where a is the angle between the rays of the sun and the longitudinal 
axis as measured with the sun camera. Assumed values of a z , a x , and a 
for various times during the calibration are given in figure 1 and 
table 1. The computations for determining a v are also indicated in 
table 1 . 

The assumed variation of p T , p^, ari T m with time is shown in 

figure 2 and table 2. Computation of temperature T* for a temperature 
recovery factor of the thermometer of 0.99 is also shown in table 2. 

The initial vertical velocity v^ is computed from the data 
assumed for the statoscope pressure recorder as shown in figure 3 . The 
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static— pressure error ep at the start of the dive is assumed to be 

known and to have the value 0.02. With the use of this value of static- 
pressure error, the statoscope data may be corrected as follows: 


Ap = Apm( 1 - e p ) 


This correction, however, is negligible for the present example. From 
a plot of Ap against t, the slope dp/dt at t = 0 is found to be 
0.04 inch of water per second. The initial vertical velocity . v-^ is 

then found from the expression 



where 


T» x = 39^-^ (table 2) 

P 1 = p mi ( 1 “ e p) ( table 2 ) 

= 76.8 (l — 0.02) =75.3 inches of water 

and 


R = 53.3 


Therefore 


= —11.2 feet per second 

Pm " p 

The computation of the static— pressure error r of the 

P 

airspeed installation is indicated in table 3* The static— pressure 
error is also plotted in figure b against indicated free— stream Mach 
number . 
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TABLE 2.- COMPUTATION OF INDICATED MACH NUMBER AND APPROXIMATE 
FREE- STREAM TEMPERATURE FOR HYPOTHETICAL CALIBRATION 
FROM PRESSURE AND TEMPERATURE DATA 

[All pressures in inches of water and temperature in °F abs?] 


t 

P«jl 

(a) 

Pm 

(a) 

P T 

Pm 

M' 

T’ 

(b) 

T m 

(a) 

T' 

0 

96.2 

76.8 

1.253 

0.577 

1.066 

420.4 

394.4 

1 

96.2 

76.8 

1.253 

.577 

1.066 

420.4 

394.4 

2 

96.2 

76.9 

1.251 

.575 

1.066 

420.4 

394.4 

3 

96.2 

76.9 

1.251 

.575 

1.066 

420.4 

394.4 

4 

96.3 

77-1 

1.249 

.573 

1.065 

420.4 

394.7 

5 

97.0 

77-3 

1.255 

• 579 

1.066 

420.7 

394.7 

6 

97-9 

77-6 

1.262 

. 586 

1.068 

421.3 

394.5 

7 

99-0 

78.1 

1.268 

• 592 

1.069 

421.9 

394.7 

8 

100.2 

78.7 

1.273 

.598 

1.071 

422.6 

394.6 

9 

101.8 

79.4 

1.282 

.607 

1.073 

423.5 

394.7 

10 

103.8 

80.2 

1.294 

.618 

1.076 

424. 5 

394; 5' 

ll 

106.1 

81.2 

1.307 

.630 

1.079 

425.5 

394.3 

12 

108.8 

82.5 

1.319 

.642 

1.082 

426.8 

394.5 

13 

111.9 

83.8 

1.335 

.656 

1.085 

428.3 

394.7 

14 

115-5 

85.4 

. 1.352 

.671 

1.089 

430.0 

394.9 

15 

119-5 ' j 

87.O 

1-374 

.689 

1.094 

431.8 

394.7 

16 

124.3 

83.9 

1.398 

.709 

1.100 

434.1 

394.6 

17 

129.5 

90.7 

1.428 

.732 

1.106 

436.8 

394.9 

18 

135-1 

92.9 

1.454 

.751 

1.112 1 

439.7 

395.4 

19 

l4i . o- 

95-0 

1.464 

• 773 

1.118 

442.9 

396.2 

20 

147.0 

97.4 

1.509 

.790 

1.124 

. 446.3 

397.1 

21 

153.3 

99-9 

1.535 

.807 

1.129 

450.0 

398.6 

22 

160.0 

102.4 

1.563 

.825 

1.135 

454.0 

4oo.o 

23 

167.2 

105-0 

1.592 

.843 

l.l4l 

458.2 

401.6 

2b 

174.6 

107.8 

1.620 

.860 

1.146 

462.9 

403-9 

25 

182.5 

110.6 

1.650 

.877 

1.152 

467-6 

405-9 

26 

190.8 

113-6 

1.680 

.894 

1.158 

472.5 

4o8.0 

27 

199.0 

116.7 

1.705 

• 907 

1.163 

476.9 

410.1 

28 

207.0 

120.0 

1.725 

.918 

1.167 

481.1 

412.3 

29 

214.7 

123-3 

1.741- 

.926 

1.170 

485.0 

414.5 

30 

222.1 

126.6 

1.754 

.933 

I.172 

488.5 

4l6.8 

31 

229-5 

130.0 

1.765 

• 939 

1.175 

491.8 • 

4l8.6 

32 

236.8 

133.5 

1.774 

.943 

1.176 

495.0 

420.9 

33 

243.7 

137-0 

1.779 

.946 

1.177 

497.8 

422.9 

3 b 

250.5 

i4o.6 

1.782 

.947 

1.178 

500.5 

424.9 


a From recorded data. See also figure 2. 
b — = 1 + 0.2KM' 2 = 1 + 0.198M' 2 where 


K = 0.99- 




TABLE 3.- CALCULATION OF FREE— STREAM STATIC PRESSURE AND STATIC-PRESSURE ERROR 
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Figure 3.- Assumed variation with time of change in static pressure at 
start of hypothetical calibration. Apjjj = Pm — p^ . 
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Figure 



w - P, deg 

(a) Due to error'of 110 miles in latitude. 



(b) Due to error of 110 miles in longitude. 


Error in attitude angle due to error in estimating position 
of airplane. e = 15 °; T = ^0°. 
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0 20 40 60 80 100 

9, deg 

(a) Due to error in longitudinal acceleration. 



(b) Due to error in normal acceleration. 


Figure 8.- Error in vertical component of acceleration due to an error 
of ±0.01g in longitudinal and normal accelerations. 
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Figure 12.- Error in vertical velocity due to error of ±0.01 inch 
of water per second in time rate of change of static pressure. 


At, sec 


Figure 13 . - Error in static pressure due to error of ±0.01 inch of 
. water per second in time rate of change of static pressure, 
h « 40,000 feet. 
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For error in p*p 


For error in T y 


Figure 14.- Error in free-stream static pressure due to ±l-percent error 
in measured temperature and ±l-percent error in total pressure. 


Figure 15.- Error in Mach number corresponding to ±1 -percent error 

in static pressure. * - 
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